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The thermal conductivity of the color-flavor locked phase of dense quark matter 
is calculated. The dominant contribution to the conductivity comes from photons 
and Nambu-Goldstone bosons associated with the breaking of baryon number, both 
of which are trapped in the quark core. Because of their very large mean free path 
the conductivity is also very large. The cooling of the quark core arises mostly 
from the heat flux across the surface of direct contact with the nuclear matter. As 
the thermal conductivity of the neighboring layer is also high, the whole interior 
of the star should be nearly isothermal. Our results imply that the cooling time 
of compact stars with color-flavor locked quark cores is similar to that of ordinary 
neutron stars. 



1. Introduction 

At sufficiently high baryon density the nucleons in nuclear matter should 
melt into quarks so that the system becomes a quark liquid, ft should 
be weakly interacting due to asymptotic freedom [1], however, it cannot 
be described as a simple Fermi liquid. This is due to the nonvanishing 
attractive interaction in the color antitriplet quark-quark channel, provided 
by one-gluon exchange, which renders a highly degenerate Fermi surface 
unstable with respect to Cooper pairing. As a result the true ground state 
of dense quark matter is, in fact, a color superconductor [2]. 

Recent phcnomenological [3] and microscopic studies [4,5,6,7] have con- 
firmed that quark matter at a sufficiently high density undergoes a phase 
transition into a color superconducting state. Phenomenological studies 
are expected to be appropriate to intermediate baryon densities, while mi- 
croscopic approaches are strictly applicable at asymptotic densities where 
perturbation theory can be used. It is remarkable that both approaches 
concur that the superconducting order parameter (which determines the 
gap A in the quark spectrum) lies between 10 and 100 MeV for baryon 
densities existing in the cores of compact stars. 
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At realistic baryon densities only the three lightest quarks can partici- 
pate in the pairing dynamics. The masses of the quarks are much smaller 
than the baryon chemical potential, thus, to a good approximation, all 
three flavors participate equally in the color condensation. The ground 
state is then the so-called color flavor locked (CFL) phase [8]. The original 
gauge symmetry SU{3)c and the global chiral symmetry SU{3)l x SU{3)r 
break down to a global diagonal "locked" SU{3)c+l+r subgroup. Because 
of the Higgs mechanism the gluons become massive and decouple from the 
infrared dynamics. The quarks also decouple because large gaps develop 
in their energy spectra. The breaking of the chiral symmetry leads to the 
appearance of an octet of pscudo-Nambu-Goldstonc (NG) bosons (tt^, tt^, 
K"^, K^, Tj). In addition an extra NG boson (p and a pseudo-NG boson 
7]' appear in the low energy spectrum as a result of the breaking of global 
baryon number symmetry and approximate U{1)a symmetry, respectively. 

The low energy action for the NG bosons in the limit of asymptotically 
large densities was derived in Refs. [9,10]. By making use of an auxiliary 
"gauge" symmetry, it was suggested in Ref. [11] that the low energy action 
of Refs. [9,10] should be modified by adding a time-like covariant derivative 
to the action of the composite field. Under a favorable choice of parameters, 
the modified action predicted kaon condensation in the CFL phase. Some 
unusual properties of such a condensate were discussed in Ref. [12]. 

While the general structure of the low energy action in the CFL phase 
can be established by symmetry arguments alone [9], the values of the pa- 
rameters in such an action can be rigorously derived only at asymptotically 
large baryon densities [10,11]. Thus, in the most interesting case of inter- 
mediate densities existing in the cores of compact stars, the details of the 
action are not well known. For the purposes of the present paper, how- 
ever, it suffices to know that there are 9 massive pseudo-NG bosons and 
one masslcss NG boson (p in the low energy spectrum. If kaons condense 
[11] an additional NG boson should appear. These NG bosons should be 
relevant for the kinetic properties of dense quark matter. 

It has been found [13] that neutrino and photon emission rates for 
the CFL phase are very small so that they would be inefficient in cooling 
the core of a neutron star. The purpose of the present investigation is 
to determine quantitatively the thermal conductivity of the CFL phase of 
dense quark matter in order to see whether it can significantly impact the 
cooling rate. We shall argue that the temperature of the CFL core, as 
well as the neighboring neutron layer which is in contact with the core, 
falls quickly due to the very high thermal condiictivities on both sides of 
the interface. In fact, to a good approximation, the interior of the star 
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is isothermal. A noticeable gradient of the temperature appears only in a 
relatively thin surface layer of the star where a finite flux of energy is carried 
outwards by photon diffusion. A more complete account of this work is to 
be found in Ref. [14]. 

2. Thermal conductivity 

A detailed understanding of the cooling mechanism of a compact star with 
a quark core is not complete without a study of thermal conductivity effects 
in the color STipc;rconducting quark core. The conductivity, as well as the 
other kinetic properties of quark matter in the CFL phase, is dominated by 
the low energy degrees of freedom. It is clear then that at all temperatures 
of interest to us, T <C A, it is crucial to consider the contributions of the NG 
bosons. In addition, there may be an equally important contribution due to 
photons; this is discussed in Sec. 4. Note that, at such small temperatures, 
the gluon and the quark quasiparticles become completely irrelevant. For 
example, a typical quark contribution to a transport coefficient would be 
exponentially suppressed by the factor exp(— A/T). 

Let us start from the general definition of the thermal conductivity 
as a characteristic of a system which is forced out of equilibrium by a 
temperature gradient. In response to such a gradient transport of heat is 
induced. Formally this is described by the following relation: 



where Ui is the heat current, and n is the heat conductivity. As is clear 
from this relation, the heat flow would persist until a state of uniform 
temperature is reached. The higher the conductivity, the shorter the time 
for this relaxation. 

In the linear response approximation, the thermal conductivity is given 
in terms of the heat current correlator by a Kubo-type formula. We derive 
the expression for the heat (energy) current carried by a single (pseudo-) 
NG boson field ^p. The corresponding Lagrangian density reads 



where the ellipsis stand for the self-interaction terms as well as interactions 
with other fields. Notice that we introduced explicitly the velocity param- 
eter V. In microscopic studies of color superconducting phases, which are 
valid at very large densities, this velocity is equal to I/a/S for all (pseudo-) 
NG bosons. It is smaller than 1 because Lorentz symmetry is broken due 
to the finite value of the quark chemical potential. By making use of the 



Ui = -KdiT, 



(1) 



L = - {doipdo^ - v^dupdup - m^ip^) + ... 



(2) 
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above Lagrangian density, we derive the following expression for the heat 
current: 

dL 



2T fT^b n 
where, in the Matsubara formalism, 



did'^f"'^ " 9^'^3o<^. (3) 

This definition leads to an expression for the heat conductivity in terms of 
the corresponding correlator [15]: 

n, 

4 f (Pk 

n ^ ' 

xS{inn,k)S{inn-m,k). (5) 

Here fi„ = 27rnT is the bosonic Matsubara frequency, and S{i^n, k) is the 
propagator of the (pseudo-) NG boson. In general, the propagator should 
have the following form: 

S{w, k) = i ^ , (6) 

where the width parameter r(w, k) is related to the inverse lifetime (as well 
as the mean free path) of the boson. In our calculation, it is very convenient 
to utilize the spectral representation of the propagator. 

Then the conductivity is expressed through the spectral function A{uj,k) 
as follows: 

By making use of the explicit form of the propagator in Eq. (6), we see that 
the spectral function of the (pseudo-) NG boson is 

= (a;2_e2_r2/4)2+^2r2 ' (9) 

where = ^/v'^k'^ + rv?. Because of the rotational symmetry of the system, 

the conductivity is characterized by a single scalar quantity n which is 
introduced as follows: Kij = nSij. The explicit expression for this scalar 
function reads 

48V2Tr^vTT^ Jo sinh^^V^ J ' ' 
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where we introduced the notation X = lo'^ — rn^ — For our purposes it 

will be sufficient to consider the conductivity in the limit of a small width, 
r — > 0. This is because the (pseudo-) NG bosons in the CFL quark matter 
are weakly interacting. Thus, wc derive the following approximate relation: 

« = ^^A^ r {c^ - m'f . (11) 



247r2«r2r7,„ sinh^^ 
At small temperature, T <^m, this result is further approximated by 

This demonstrates clearly that the contributions of heavy pseudo-NG 
bosons to the thermal conductivity are suppressed. The largest contribution 
comes from the massless NG boson for which the thermal conductivity is 

Stt'^vT^ Jq sinh^x 4:5vT^ 

In order to calculate the width is required, or equivalently the mean 
free path £^ since = v/T^, where v is the average thermal velocity of the 
NG bosons. This will be discussed in the next section. 



'^'P = ^^23^ / -72- = 7^ ■ (13) 



3. Mean free path of the NG boson 

As we have remarked, the contribution of massive pseudo-NG bosons to 
the thermal conductivity is suppressed. In the CFL phase of quark matter, 
however, there is one truly massless NG boson ^ which should therefore give 
the dominant contribution to the heat conductivity. The interactions of 
with the CFL matter leads to a finite value for its mean free path. Since this 
boson is a composite particle there is always a non-zero probability at finite 
temperature for its decay into a pair of quark quasiparticles. It is natural 
to expect that such a process is strongly suppressed at small temperatures, 
T <C A. This is confirmed by a direct microscopic calculation in the region 
of asymptotic densities which yields a decay width [16]: 

5V2nvk ( /3A\ 

^^--^^^ " 4(21-81n2) ^"P 2 T J " ^^^^ 

If this were the only contribution, then the order of magnitude of the mean 
free path of the NG boson would be 




6 



This grows exponentially with decreasing temperature. For example, if 
A ~ 50 MeV and T < 1.5 MeV, the mean free path is 30 km or more. This 
scale is a few times larger than the typical size of a compact star. 

The decay channel of the NG bosons into quarks is not the only contri- 
bution to the mean free path. They can also scatter on one another. The 
corresponding amplitude is of order fc^Z/U^ [17] which gives a cross section 
of a^ij, ~ T^/n^, yielding the following contribution to the width: 

^4,4, = va^^Urf, ~ — , (16) 

where n<^ is the equilibrium number density of the NG bosons [14]. At 
small temperatures the scattering contribution in Eq. (16) dominates the 
width. This leads to a mean free path 

^,0-^«8xlOsi|°iLkm. (17) 

Here we defined the following dimensionless quantities: /X500 = 
/^/(500 MeV) and TmsV = T/{1 MeV). Both (.4,^ and ^^^^^ depend very 
strongly on temperature, however the salient point is that they are both 
larger than the size of a compact star for temperatures TmcV of order 1. 

We define T to be the temperature at which the massive NG bosons 
decouple from the system. This is determined by the mass of the light- 
est pseudo-NG boson for which it is not presently possible to give a reli- 
able value. Different model calculations [10,11,18] produce different values 
which can range as low as 10 MeV. Thus, conservatively, we choose T ^ 1 
MeV. Then the mean free path of the NG boson is comparable to or even 
larger than the size of a star for essentially all temperatures T < T. It is 
also important to note that the mean free path is very sensitive to temper- 
ature changes. In particular, at temperatures just a few times higher than 
T the value of I may already become much smaller than the star size. This 
suggests that, during the first few seconds after the supernova explosion 
when the temperatures remain considerably higher than T, a noticeable 
temperature gradient may exist in the quark core. This should relax very 
quickly because of the combined effect of cooling (which is very efficient 
at T ^ T) and diffusion. After that almost the whole interior of the star 
would become isothermal. 

Before concluding this section, we point out that the geometrical size 
of the quark phase limits the mean free path of the NG boson since the 
scattering with the boundary should also be taken into account. It is clear 
from simple geometry that i Rq, where Rq is the radius of the quark 
core. 
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4. Photon contributions 

Now, let us discuss the role of photons in the CFL quark core. It was argued 
in Ref. [13] that the mean free path of photons is larger than the typical 
size of a compact star at all temperatures T < T. One might conclude 
therefore that all photons would leave the stellar core shortly after the 
core becomes transparent. If this were so the photons would be able to 
contribute neither to the thermodynamic nor to the kinetic properties of the 
quark core. However the neighboring neutron matter has very good metallic 
properties due to the presence of a considerable number of electrons. As is 
known from plasma physics, low frequency electromagnetic waves cannot 
propagate inside a plasma. Moreover, an incoming electromagnetic wave is 
reflected from the surface of such a plasma [19]. In particular, if Op is the 
value of the plasma frequency of the nuclear matter, then all photons with 
frequencies w < Vlp are reflected from the boundary. This eff'ect is similar 
to the well known reflection of radio waves from the Earth's ionosphere. 

The plasma frequency is known to be proportional to the square root of 
the density of charge carriers and inversely proportional to the square root 
of their mass. It is clear therefore that the electrons, rather than the more 
massive protons, will lead to the largest value of the plasma frequency in 
nuclear matter. Our estimate for the value of this frequency is 

^ J^f^^ 4.7 ^ ,q2 [K MeV, (18) 
y niemp y po 

where the electron density Ug = Yep/nip is given in terms of the nuclear 
matter density p and the proton mass m^. Also m,e denotes the electron 
mass, Yg ~ 0.1 is the number of electrons per baryon, and po « 2.8 x 
10^^ g cm~^ is equilibrium nuclear matter density. 

Since Qp is more than 100 MeV, essentially all thermally populated 
electromagnetic waves at T < T will be reflected back into the core region. 
In a way the boundary of the core looks like a good quality mirror with 
some leakage which will allow a thermal photon distribution to build up 
and stay. Thus photons will be trapped in such a core surrounded by a 
nuclear layer. Notice that the transparency of the core is reached only after 
the temperature drops substantially below T ~ 2me, i.e., when the density 
of thermally excited electron-positron pairs becomes very small. 

Now, since photons are massless they also give a sizable contribution to 
the thermal conductivity of the CFL phase. The corresponding contribution 
Kj will be similar to the contribution of massless NG bosons in Eq. (13). 
Since the photons move at approximately the speed of light [20] (v ~ 1 at 
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the densities of interest) and they have two polarization states, we obtain 

Since the thermal conductivity is additive the total conductivity of dense 
quark matter in the CFL phase is given by the sum of the two contributions: 

I^CFL = + ^ ——T^Ro, (20) 
y 

where for both a photon and a NG boson the mean free path i ^ Rq. This 
yields the value 

KcFL - 1.2 X 1032T3i^yi?o,k„ erg cm-isec-iK-i, (21) 

where i?o,km is the quark core radius measured in kilometers. The value of 
KcFL is many orders of magnitude larger than the thermal conductivity of 
regular nuclear matter in a neutron star [21]. 



5. Stellar cooling 

In discussing the cooling mechanism for a compact star wc have to make 
some general assumptions about the structure of the star. We accept with- 
out proof that a quark core exists at the center of the star. This core stays 
in direct contact with the neighboring nuclear matter. From this nuclear 
layer outwards the structure of the star is essentially the same as an or- 
dinary neutron star. The radius of the core is denoted by Rq, while the 
radius of the whole star is denoted by R. 

A detailed analysis of the interface between the quark core and the 
nuclear matter was made in Ref. [22]. A similar analysis might also be 
very useful for understanding the mechanism of heat transfer from one 
phase to the other. Here we assume that direct contact between the phases 
is possible, and that the temperature is slowly varying across the interface. 

Now, let us consider the physics that governs stellar cooling. We start 
from the moment when the star is formed in a supernova explosion. Imme- 
diately after the explosion many high-energy neutrinos arc trapped inside 
the star. After about 10 to 15 seconds most of them escape from the star 
by diffusion. The presence of the CFL quark core could slightly modify the 
rate of such diffusion [23,24,25]. By the end of the deleptonization process, 
the temperature of the star will have risen to a few tens of McV. Then, 
the star cools down relatively quickly to about T by the efficient process of 
neutrino emission. It is unlikely that the quark core would greatly affect 
the time scale for this initial cooling stage. An ordinary neutron star would 
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continue to cool by neutrino emission for quite a long time even after that 
[26]. Here we discuss how the presence of the CFL quark core affects the 
cooling process of the star after the temperature drops below T. 

Our result for the mean free path of the NG boson demonstrates that 
the heat conductivity of dense quark matter in the CFL phase is very high. 
For example, a temperature gradient of 1 MeV across a core of 1 km in size 
is washed away by heat conduction in a very short time interval of order 
iig y^^/v(.(T) ~ 6 X 10"** sec. In deriving this estimate, we took into account 
the fact that the specific heat and the heat conductivity in the CFL phase 
are dominated by photons and massless NG bosons and that ^ ~ i?o- In 
addition, we used the classical relation k = vc^,C./'.i, where is the specific 
heat; this can be shown to hold in the present context [14]. Since heat 
conduction removes a temperature gradient in such a short time interval, 
it is clear that, to a good approximation, the quark core is isothermal at 
all temperatures T <T. 

The heat conductivity of the neighboring nuclear matter is also known 
to be very high because of the large contribution from degenerate electrons 
which have a very long mean free path. It is clear, then, that both the 
quark and the nuclear layers should be nearly isothermal with equal values 
of the temperature. When one of the layers cools down by any mechanism, 
the temperature of the other will adjust almost immediately due to the very 
efficient heat transfer on both sides of the interface. 

Now consider the order of magnitude of the cooling time for a star with a 
CFL quark core. One of the most important components of the calculation 
of the cooling time is the thermal energy of the star which is the amount of 
energy that is lost in cooling. There are contributions to the total thermal 
energy from both the quark and the nuclear parts of the star. The dominant 
amount of thermal energy in the CFL quark matter is stored in photons 
and massless NG bosons. Numerically, its value is [14] 

Ecfl{T) ^ 2.1 X lO^^^o'.kmr^icV erg. (22) 

The thermal energy of the outer nuclear layer is provided mostly by degen- 
erate neutrons. The corresponding numerical estimate is [27]: 

E.m{T) ^ 8.1 X 10''^^ (^) T^ev erg, (23) 

where M is the mass of the star, Mq is the mass of the quark core and Mq 
is the mass of the Sun. It is crucial to note that the thermal energy of the 
quark core is negligible in comparison to that of the nuclear layer. 

The second important component that determines stellar cooling is the 
luminosity which describes the rate of energy loss due to neutrino and 
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photon emission. Typically, the neutrino luminosity dominates the cooling 
of young stars when the temperatures are still higher than about 10 keV 
and after that the photon diffusion mechanism starts to dominate. Photon 
and neutrino emission from the CFL quark phase is strongly suppressed at 
low temperatures [13]. The neighboring nuclear layer, on the other hand, 
emits neutrinos quite efficiently. As a result, it cools relatively fast in the 
same way as an ordinary neutron star. The nuclear layer should be able to 
emit not only its own thermal energy, but also that of the quark core which 
constantly arrives by the very efficient heat transfer process. The analysis 
of the cooling mechanism is greatly simplified by the fact that the thermal 
energy of the quark core is negligible compared to the energy stored in the 
nuclear matter. By taking this into account, we conclude that the cooling 
time of a star with a quark core is essentially the same as for an ordinary 
neutron star provided that the nuclear layer is not extremely thin. 

6. Conclusions 

Our analysis shows that the thermal conductivity of CFL color supercon- 
ducting dense quark matter is very high for typical values of the temper- 
ature found in a newborn compact star. This is a direct consequence of 
the existence of the photon and the massless NG boson (p whose mean free 
paths are very large. Note that the photons are trapped in the core because 
of reflection by the electron plasma in the neighboring nuclear matter. The 
NG bosons are also confined to the core since they can only exist in the 
CFL phase. 

It is appropriate to mention that the (pseudo-) NG bosons and photons 
should also dominate other kinetic properties of dense quark matter in 
the CFL phase. For example, the shear viscosity should be mostly due to 
photons and the same massless NG bosons associated with the breaking of 
baryon number. The electrical conductivity, on the other hand, would be 
mostly due to the lightest charged pseudo-NG boson, i.e., the K^. Thus, 
in the limit of small temperatures, T ^ 0, the electrical conductivity will 
be suppressed by a factor exp(— mx+/r). 

Since the neutrino cmissivity of the CFL core is strongly suppressed, the 
heat is transferred to the outer nuclear layer only through direct surface 
contact. While both the core and the outer layer contribute to the heat 
capacity of the star, it is only the outer layer which is capable of emitting 
this heat energy efficiently in the form of neutrinos. The combination of 
these two factors tends to extend the cooling time of a star. However, 
because of the very small thermal energy of the quark core, the time scale 
for cooling could be noticeably longer than that for an ordinary neutron 
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star only if the outer nuclear layer was very thin. (Note that, while little is 
known about the properties of thin boundary layers outside CFL matter, 
it is possible that photon emission from this layer might quickly drain the 
relatively small amount of CFL thermal energy.) Thus it appears that the 
cooling of stars with not too large CFL quark cores will differ little from 
the cooling of typical neutron stars. A similar conclusion has been reached 
for stars with regular, non-CFL quark interiors [28]. 

In passing it is interesting to speculate about the possibility that a bare 
CFL quark star made entirely of dense quark matter could exist. If it 
were possible, it would look like a transparent dielectric [29]. Our present 
study suggests such a star would also have very unusual thermal properties. 
Indeed, if the star has a finite temperature T < T after it was created, 
almost all of its thermal energy would be stored in the NG bosons. Notice 
that all the photons would leave the star very soon after transparency set 
in because the star is assumed to have no nuclear matter layer. The local 
interaction as well as the self-interaction of the NG bosons is very weak so 
that we argued in Sec. 3 that their mean free path would be limited only 
by the geometrical size of the star. This suggests that, since photon and 
neutrino emission inside the CFL phase is strongly suppressed [13], the only 
potential source of energy loss in the bare CFL star would be the interaction 
of the NG bosons at the stellar boundary and photon emission there. It is 
likc^ly. therefore, that such stars might be very dim and might even be good 
candidates for some of the baryonic dark matter in the Universe. 
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